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Abstract 

A new upper limit on the ampHtude of primordial magnetic field (PMF) is derived by a com- 
parison between a calculation of elemental abundances in big bang nucleosynthesis (BBN) model 
and the latest observational constraints on the abundances. Updated nuclear reaction rates are 
adopted in the calculation. Effects of PMF on the abundances are consistently taken into account 
in the numerical calculation with the precise formulation of changes in physical variables. We find 
that abundances of ^He and ^Li increase while that of ^Li decreases when the PMF amplitude 
increases, in the case of the baryon-to-photon ratio determined from the measurement of cosmic 
microwave background radiation. We derive a constraint on the present amplitude of PMF, i.e., 
-6(0) < 1.5 /xG [corresponding to the amplitude less than 2.0 x 10^^ G at BBN temperature of 
T = 10^ K] based on the rigorous calculation. 

PACS numbers: 26.35. +c, 98.62.En, 98.80.Es, 98.80.Ft 



* |kusakabe@icrr . u-tokyo .ac.jp 



1 



I. INTRODUCTION 



Primordial nucleosynthesis, or big bang nucleosynthesis (BBN), has been assumed [l| to 
occur through complicated nonequilibrium processes. It involves many reactions including 
radiative neutron capture reactions l|-3| and weak interactions converting protons and neu- 
trons to each others j4| as well as relativistic quantum statistics jsl- In BBN, only D, ^He, 
^He and ^Li can be produced in significant amounts jo], and yields of heavier elements are 
generally expected to be small 0, IsJ. The BBN model predicts the relic of a dense hot radi- 
ation j^, 9] to be observed as cosmic microwave background radiation (CMBR) today 10]. 

The BBN has been studied over a long period, and the theory is now precisely struc- 
tured (e.g., jol. [llijssll ) . The simplest, standard BBN (SBBN) model is characterized by one 
parameter, i.e, baryon-to-photon number ratio t] with the fixed number of light neutrino 
species of = 3. The t] value is constrained precisely with data of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 39l-l4l|. The SBBN model prediction of light element 
abundances for the WMAP r] value is rather consistent with primordial abundances inferred 
from observations. There is, however, a discrepancy between the predicted and observed 
primordial abundances of ^Li. The SBBN predicts a ^Li abundance which is a factor of 



2.4 — 4.3 times higher 



42] than the observationally deduced abundance. Possible solutions 



to this discrepancy have been proposed (e.g. 43| and references therein). 



O'Connel and Matese 4^ have estimated the neutron /3-decay rate in the presence of a 
strong magnetic field, and suggested that an increase in the rate due to primordial magnetic 
fields (PMFs) decreases '^He abundance. Greenstein 45| subsequently suggested that the 
energy of PMFs enhances the expansion rate of the universe, and it tends to increase the 
^He abundance rather than decrease it as suggested in Ref. 4J]. Matese & O'Connel 46 ] 
then performed a detailed investigation on the PMF effects on BBN, and concluded that the 
effect through the expansion rate is predominant over that through rates of weak reactions. 



47H52], and have the common 



Three groups have investigated effects of PMF on BBN 
conclusion that the effect through the cosmic expansion rate contributed from an enhanced 
energy density 45| is the most important 49l-l5l|. The effect of energy density of PMF 
can be considered in analogy with that of an effective neutrino number in BBN epoch j49l |. 
Grasso & Rubinstein 50| have additionally shown that a change in the quantum statistics 



of electron and positron by the PMF affects BBN. Constraints on PMF depend on other 



parameters than the amphtude of PMF. Suh & Mathews [53| have studied sensitivity of 



hmits on PMF to the neutrino degeneracy. See section 3 in Ref. 
topic. 



54j for a review of this 



The latest BBN constraint on the magnetic field 51 



has been vary old. It was 
based on assumptions of an old baryon-to-photon ratio rj = 2.8 x 10~^° and the upper limit 
on ^He mass fraction of Yp < 0.245. The values are updated to be r/ = 6.2 x 10"^*^ for 



ACDM+SZ+lens model 



41| . and Y = 0.2561 ± 0.0108 55|. In this paper we perform a 



network calculation of BBN taking account of effects of PMFs, and show a latest constraint 
on PMF through effects on elemental abundances. In addition, formulae necessary for 
precise numerical calculations are provided. This study improves the following points over 
previous works: 1) updates on nuclear reaction rates, the neutron lifetime, and observational 
constraints on primordial abundances, 2) a precise treatment of electron chemical potential 
in abundance calculation and an estimation for initial value of electron chemical potential, 
3) a precise calculation of temperature evolution as a function of time or cosmic scale factor, 
and 4) a caution that an effect of magnetic field on nuclear reaction rates is weak. 

In Sec. iniwe describe the model of SBBN code with a recent update on nuclear reaction 
rates and the neutron lifetime, and also how to include magnetic fields effects on BBN in 
precise numerical studies. In Sec. Illll we show results of calculations of BBN in the presence 
of variable amplitudes of PMF. In Sec. IIVI we discuss constraints on PMFs. In Sec. |V] 
we summarize this study. In Appendix El we describe formulae necessary for BBN network 
calculations including effects of PMFs. In Appendix |B] an effect of PMFs on nuclear reaction 
rates is studied, and it is shown to be negligible. 



II. MODEL 



A. standard BBN 



We use a BBN code Il9l . 1561 1 for reaction network calculations. The Sarkar's correction is 
adopted for ^He abundance 



23| . Rates and their uncertainties of reactions for light nuclei 



{A < 10) are updated with recommendations of the JINA REACLIB Database VI. (571. We 
derive 95 % confidence regions of elemental abundances assuming uncertainties in rates of the 



12 important reactions 



19| . The rates are assumed to be given by the Gaussian distribution. 



TABLE I. Adopted reaction rates 





reaction 


reference 


1 


n(,e"z7g)iH 


[58] and [59] 


12 


iH(n,7)2H 


[60] 


16 


3He(n,p)3H 


[6^ 


17 


'^Be(n,p)^Li 


[61J 


20 


2H(p,7)3He 


[61J 


24 


^Li(p,Q)4He 


[61J 


26 


3H(a,7)^Li 


m 


27 


3He(a,7)7Be 


[62] 


28 


2H(d,n)3He 


[61J 


29 


2H((i,p)3H 


[61J 


30 


3H(d,n)4He 


[61J 


31 




m 



Reaction number in tlic Kawano's code 



and 1000 runs are performed for each eta value. The reactions and the references for adopted 
rates are hsted in Table HI 

We adopt two values of neutron lifetime. One is 878.5 ± 0.7stat ± O.Sgys s from Ref. [gS] 
based on improvements jssj in the measurement. This relatively short lifetime better satisfies 

and it can improve the 



64 



65 



Another 



59j. As of 



the unitarity test of the Cabibbo-Kobayashi-Maskawa matrix [58] 
agreement between observed primordial abundances and BBN predictions 
value is 885.7 ± 0.8 s from the old recommendation by the Particle Data Group 
April 2012, the Particle Data Group presents a new average neutron lifetime of 881.5 ± 1.5 s 
which is sandwiched between the adopted lifetimes. jfiS]. 

We adopt constraints on primordial abundances as follows: 

A deuterium abundance in a damped Lyman alpha system of QSO SDSS J1419+0829 
was measured precisely than any other QSO absorption systems (67|. We adopt both of a 
mean value of ten QSO absorption line systems including J1419+0829, and the abundance of 
J1419+0829 itself, i.e., log(D/H)=-4.58 ± 0.02 and log(D/H)=-4.596 ± 0.009, respectively. 



We take 2a uncertainties, i.e., 

2.40 X 10"^ < D/H < 2.88 X 10"^ (mean), 

2.43 X 10"^ < D/H < 2.64 X 10"^ (best). (1) 

■^He abundances are measured in Galactic HII regions through the 8.665 GHz hyperfine 
transition of '^He"'", i.e., ^He/H=(1.9 ± 0.6) x 10~^ 68|]. Although the constraint is rather 
weak considering its uncertainty, we take a 2cr upper limit from abundances in Galactic HII 
region, i.e., 

^He/H < 3.1 X 10"^ (2) 
For the primordial helium abundance we adopt two different constraints, i.e, Y = 0.2565± 



0.0051 



69j and F = 0.2561 ± 0.0108 



both from observations of metal-poor extragalactic 
HII regions. We take 2a limits of 

0.2463 <Y < 0.2667 (ITIO), 

0.2345 <Y < 0.2777 (AOSIO). (3) 

As a guide, observed lithium abundances follow although they are not used as constraints. 

Primordial ^Li abundance is inferred from spectroscopic observations of metal-poor halo 
stars (MPHSs). We adopt log(^Li/H)= -12 + (2.199 ± 0.086) (95% confidence limits) 
derived in a 3D nonlocal thermal equilibrium model 70|]. This estimation corresponds to 
the 2a range of 

1.06 X 10"^° < ^Li/H < 2.35 x 10"^°. (4) 

le stars. The most 



Tlj. We use the 2a 



Observations of MPHSs suggest a presence of ^Li nuclei in some of t 
probable detection of ^Li for G020-024 indicates ''Li/^Li=0.052 ± 0.017 
upper limit and log(^Li/H) = —12 + 2.18 for the same star 72|, and derive 

^Li/H < 1.3 X 10^". (5) 

Figure [1] shows abundances of ^He [Yp] mass fraction), D, "^He, ^Li and ^Li (A/H; by 
number relative to H) as a function of the baryon-to-photon ratio t] or the baryon energy 
density fi^/i^ of the universe. The solid and dashed curves are the results for neutron 



lifetimes of 878.5 ± 0.8 s 58| and 885.7 ± 0.8 s 59|], respectively. Thin solid curves show 95 
% ranges determined from uncertainties in nuclear reaction rates. The boxes represent the 
adopted abundance constraints. The vertical stripe represents the 2 a Vtsh'^ limits provided 



by WMAP 



41]. This corresponds to fis/i^ = 0.02258l[]:[]|]}i^ or r/ = (6.225lgJJ^) x 10" 
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FIG. 1. Abundances of '^He (mass fraction), D, "^He, ^Li and ^Li (by number relative to H) as a 



function of the baryon-to-photon ratio r] or the baryon energy 



density f^s/i^. The so. 



581] and 885.7 ±0.8 s 



curves are the results for neutron lifetimes of 878.5 it 0.8 s 

Thin solid curves show 95 % ranges determined from uncertainties in nuclear reaction rates. 

691 for ^He, 



id and dashed 
59|, respectively. 

he 



boxes represent the adopted observational abundances from Refs. 
for 3He, Q for ^Li, and 71 



limits provided by WMAP 



551, 



67|] for D, 



68| 



72l | for ^Li, respectively. The vertical stripe represents the 2 a ri^/i^ 



ill], i.e., Qsh^ = 0.02258 



+0.00114 
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B. effects of magnetic field 

We include effects of PMF through the magnetic energy density (Appendix [1]) , thermo- 
dynamic variables of electron and positron and their time evolutions (Appendixes [21 |3]). 
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Estimations for initial values of electron chemical potential are changed from those in the 
case of no magnetic field (Appendix H]) . Equations to solve are similar to those in Ref. 49|. 
However, equations which are solved in the consistent numerical calculation (see Appendix 
R|) are more complicated. 

Final r] values are different for different initial B values with a fixed initial t] value. The 
final r] value, instead of the initial value, should then be fixed to the value determined from 



WMAP 



4l| as pointed out but not done in deriving the limit on B [eq. (4)] in Ref. 50 1. 



In this study we fixed the final rj value. The effect o 



the magnetic field on weak reaction 



rates has been long since found to be negligible j49|, |50|. It is, therefore, not included in this 
calculation. 



III. RESULT 



Figure [2] shows light element abundances as a function of the magnetic field amplitude in 
units of the critical value, i.e., 7 = B/Bq at temperature T = 10^ K, or absolute value at the 
present epoch of 2; = 0, i.e., 5(0). Bq = 4.41 x 10^'^ G is the critical magnetic field (above 



which quantized magnetic levels appear 5J]). The two parameters are related by 7(T = 10^ 
K) = 3.05[-B(0)/mG]. The solid and dashed curves correspond to results for two different 
neutron lifetimes, and the boxes represent adopted abundance constraints (see Sec. Ill A|) . 

Figure El shows time evolutions of light element abundances as a function of the temper- 
ature Tg = T/(10^ K). Solid lines correspond to the case of a magnetic field of B(z = 0) = 
5 fiG, while dashed lines correspond to SBBN. 

The primordial abundance of ^He increases when the amplitude of PMF increases. The 
cosmic expansion rate is larger because of the energy of the PMF, so that the neutron 
abundance after the earlier freeze-out of weak reactions is higher. The time interval between 
the freeze-out and the ^He synthesis is also shorter because of faster cosmic expansion. 
Neutron abundances are larger than in SBBN for the above two reasons. Those neutrons 
are processed to form ^He nuclei. This is the reason of the trend of ^He abundance vs. B 
value. 

Because of the earlier freeze-out of the reaction ^H(n, 7)^H at the ^He synthesis, the 
relic neutron abundance is higher than in SBBN. This higher neutron abundance affects 
abundances of other light nuclei complicatedly. The abundance of D which is produced 
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FIG. 2. Abundances of '^He (mass fraction), D, "^He, ^Li and ^Li (by number relative to H) as a 
function of the magnetic field amplitude in units of the critical value, i.e., 7 = B/Bq, or absolute 



value at the present epoc! 
lifetimes of 878.5 ± 0.8 s 



1 of z = 0. The solid and dashed curves are the results for neutron 



5i and 885.7 ± 0.8 s 



59|, respectively. The boxes represent adopted 



abundance constraints which are the same as in Fig. [TJ The final value of baryon-to-photon ratio 

n 

is fixed to be ry = 6.2 x 10"^° [41,]. 



via ^H(n, 7)^H is somewhat higher when B is higher. is produced via '^B.{d,p)^B. and 
destroyed via '^H((i, n)^He. The enhanced D abundance simply leads to a higher abun- 
dance by a higher production rate. '^He, on the other hand, is produced via ^H((i, n)^He and 
destroyed via '^He(n,p)^H. The somewhat higher D abundance leads to higher production 
rate while the higher neutron abundance leads to higher destruction rate. Resultingly, the 
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FIG. 3. Abundances of H and ^He, i.e., X and Y, respectively, (mass fraction), and other nuclides 
(by number relative to H) as a function of the temperature Tg = T/(10^ K). Solid lines correspond 
to the case of a magnetic field of B{z = 0) = 5 fxG, while dashed lines correspond to standard 
big bang nucleosynthesis model. In both cases, final values of baryon-to-photon ratios are rj = 



6.2 X 10-1° 



4l|. 



^He abundance is slightly higher than in SBBN. 

^Li is produced via ^He((i, 7)^Li and destroyed via ^Li(p, a)^He, and ''Li is produced via 
^He(t, 7)^Li and destroyed via ''Li(p, a)^He. The abundances of both nuclides are then higher 
since those of D and T are higher at higher B values. ^Be is produced via ^He(^He, 7)^Be and 
destroyed via ^Be(r2,p)^Li. Slightly higher abundance of ^He and rather higher abundance 
of neutron results in net reduction of ''Be abundance with respect to that of SBBN. 

The final ^Li abundance is the sum of those for ^Li and ''Be. ^Be nuclei convert to ^Li 
nuclei by an electron capture process. Since the abundance of ^Be is larger than that of ^Li 
in SBBN, an existence of PMF reduces the final ^Li abundance. 

Shapes of curves in Fig. [2] are explained as above. Increases in abundances of ^He and 



(D+^He) have been obtained in the previous investigation 50||, and are consistent with our 
result. 

The following constraints are derived from Fig. O 

B{0) < 2.0 i^iG or jiTg = 1) < 6.0 x lO'^ (^He; ITIO), 

5(0) < 2.5 /iG or 7(T9 = 1) < 7.7 x 10-^ (^He; AOSIO), (6) 



B{0) < 1.5 jjG or 7(T9 = 1) < 4.6 x lO'^ (D; mean), 
B{0) < 0.78 /iG or 7(T9 = 1) < 2.4 x lO^^ (D; best), 



(7) 



If one conservatively takes the constraint on ^He abundance by AOSIO, the observation 
of D abundance provides the strongest upper limits on PMF. The conservative upper limit 
from the mean value of QSO D/H ratio [Glis B{0) < 1.5 jjG [-f{Tg = 1) < 4.6 x 10"^], 
while that from the best D/H measurement [671 is B{0) < 0.78 fiG [-fiTg = 1) < 2.4 x lO'^]. 
The latter limit is nearly identical to the previous estimation (corresponding to '~f(Tg = 1) < 



2.3 X 10 which is read from eq. (4) in Ref. 50|), while 



the former by a factor of two. Previous constraints 491-151 



;he former is less stringent than 



5J] have been derived neglecting 



changes in evolution of baryon-to-photon ratio rj. In our work, this effect is consistently 
taken into account, and the final rj value is fixed to the WMAP estimation. The present 
result is, therefore, most precise. Other improvements are updates of nuclear reaction rates, 
observational constraints on primordial abundances, and baryon-to-photon ratio. 



IV. DISCUSSION 



The constraint derived in this study is related to the local field amplitude B contributed 



from all wavelengths, and is not for that measured at any particular scale 49|]. The present 



amplitude of cosmological averaged field, i.e, -Brms, (or the energy density of magnetic field) 



is defined 



49| by 



B{z)Ls = ^ I d^x\B{x,z)\\ (8) 

where z is the redshift, Vh is the Hubble volume, and x is the position vector. The conser- 
vative constraint is then i?rms(0) < 1.5 /iG. 

Magnetic fields on some scales depend on the spatial structure of field. The root mean 



square (rms) amplitude on scale L 



54 is 

T \P 



B{L,z\^, = B{<d\^,{l + zY [-^J , (9) 

where Lq is the comoving coherence length, and p is a parameter determined from statistical 
properties of the magnetic field js^ . 

10 



The present constraint can be compared with those from other measurements summarized 
m Refs HQ- 

We note that the new constraint [Eq. ([7])] can be the strongest for small 
correlation scales of Lq ^ 10"^ pc. Direct constraints on magnetic field strength at smallest 
scales are derived from observations of Zeeman effect of HI, OH and CN in molecular clouds 
and HI diffuse clouds 75|. The smallest upper limit on the radial component of magnetic 



field is B» = 0.0 ± 0.9 uG for an HI cloud seen in absorption against radio source 3C 348 



(a usable data in Ref. 76|). Heiles and Troland used data of Zeeman-splitting of the 21 cm 
line 76|, and estimated a median total field strength B = 6.0 ± 1.8 /iG for HI clouds with 



scales of typically (9(0.01 - 10 pc), taking account of probability distribution function of 
total field strength B and a random orientation of fields with respect to the line of sight [t^]. 

The constraint on PMF from BBN studies cannot be directly compared with those from 
CMBR studies (e.g. ItsHsiI) , i.e., B{1 Mpc, 0)rms = 0.85 ± 1.25 nC [Slj, since the CMBR 
limits are imposed on magnetic fields on scales larger than the horizon in the BBN epoch. 54 1. 
For example, when we adopt Lq = 100 pc (the comoving Hubble horizon in the BBN epoch) 
and p = 3/2 (which is derived in the assumption that a field vector performs a random walk 



in three dimensional space by steps of the physical scale Lq 



82|), Eq. (ED leads to 



B{1 Mpc,0)rms < 1-5 pG. 



(10) 



Comparisons between constraints for different coherent lengths thus generally depend 



on statistical properties of magnetic field. See Ref. 



for a recent review for creation 



mechanisms of extragalactic magnetic fields and their problems. 



V. CONCLUSIONS 

A new upper limit on the amplitude of primordial magnetic field (PMF) is derived by 
a comparison between a numerical calculation of elemental abundances in big bang nucle- 
osynthesis and the latest constraints on abundances inferred from observations. The newest 
nuclear reaction rates are adopted (Sec. |Tl])- In addition, effects of PMF on the abundances 
are consistently taken into account in the numerical calculation with a formulation of phys- 
ical variables in a magnetic field (Appendix IXj) . 

We find that the existence of PMF increases abundances of ^He and ^Li, and decreases 
that of ^Li in the calculation for the baryon-to-photon ratio determined from the measure- 

11 



ment of cosmic microwave background radiation with the Wilkinson Microwave Anisotropy 
Probe. As a result of the rigorous calculation, we derive a constraint on the present ampli- 
tude of PMF, i.e., -B(O) < 1.5 /iG [corresponding to the amplitude less than 4.6 x 10~^ times 
the critical magnetic field strength for electron at temperature T = 10^ K]. 



Appendix A: Formulae for effects of magnetic field on nucleosynthesis 

1. energy density 

The energy density of magnetic field is 

pj, = ^ = (Al) 

where B is the amplitude of magnetic field, and 7 = B/Bq is the B value in units 
of critical magnetic field, i.e., Bq = ml/e = 4.41 x 10^^ G with e the electric charge, 
and rrie the electron mass. This energy contributes to the total energy contents of the 
universe related to the Hubble expansion rate. In this study, it is assumed that the 
primordial magnetic field (PMF) just attenuates by the cosmic expansion. 

2. thermodynamic variables of electron in a magnetic field 



The number density, energy density and the pressure 8J| of the electron and positron 



are given 



54| . respectively, as 

eB °° r°° 
ne{B) = -—- ^(2 - 5„o) j ^fFD{Te,E^) dp,, (A2) 



eB °° /■°° 
Pe{B) = — — 5: (2 - 5„o) / En fMTe, En) dp,, (A3) 

l^TTj J-00 

Pe{B) = 7^ E(2 - 4o) fMTe.En) dp,, (A4) 



where 



l + exp[(E„(pJ =FyU) lTe\ 
is the Fermi-Dirac distribution function at electron temperature Te, and En = [pi + 
eB{2n + l + s) +mg]^/^ is the energy of electron in the presence of a uniform field which 
is much smaller than the critical strength Bq n = 0, 1, 00 and s = ±1 are 
the principal and magnetic quantum numbers of the Landau level, respectively, and p 

12 



is the chemical potential of electron. It has been assumed that the direction of field is 
the z-axis. 



The above quantities can be rewritten in the form of 



m^7 °° 



ns=0 



27r2 



ns=0 "^^ J- -r e 



(A6) 
(AT) 



dk 



ns=0 



3e 1 + e"'=T'^'= 



(A8) 



where k = pz/me, e = {k"^ + 1 + 2777.3)"^'^^, Ze = me/Te, and </)e = /i/Tg were defined. 
Using the Euler-McLaurin formula 



the number density and the energy density of 



electron and positron are given by 



7^3 I ^00 

e 



dp 



drj 



+ 



7^ /me^ ^ 



24 V Te / Jo + ^2/2^2 I + cosh(yr/2 + ^2/7^2 ^ , 

2\ 4] 



7 



(A9) 



2^4 



p^Jp'^ + ml/T^ dp 



1 + exp(w'p2 + ml/T^ T <Pe 



1 + cosh(w?72 + ml/T^ =F 



-0 



1 



T2 



1 + exp(y^?72 + ml/T^ =f 0e 

(AlO) 



These equations are the same as those derived in Ref. 49| except that ours are general- 
ized versions including the electron chemical potential. In order to follow in numerical 
calculations precisely the electron chemical potential, which becomes large at late time 
of BBN, it is kept in our formulation. Adopting the Euler-McLaurin formula to the 
pressure of electron and positron, one can obtain the equation, i.e., 

dp 



P = — 

' 37r2 



Jp^ + ml/T^ l + expUp^ + ml/T^T(l>e) 



13 



r/2 + ml/T^ I + cosh( + ml/T^ T <t>e) 



2/. + ml/T^ 



+0 



(All) 



Eqs. (IA9p . ( 1A10I) and (lAlip reproduce values for the case of no magnetic field when 
i? = is input. 

Time evolutions of following three variables as perturbations induced by a magnetic 
field are calculated. The first is related to the asymmetry in number abundances of 
electron and positron, i.e, 



TT 



he 



1 /■°° 
48 io 



where h is the Planck's constant, c is the light speed, and 



2. A - 



fniv) =1 Z 



9 



1 + cosh(6' 



1 + cosh(6' 



(A12) 



(A13) 



and the parameter 6{ri) =\l'rf + m'l/T'^ was defined. Partial derivatives of this function 
with respect to Tg = Tg/ (10^ K), the neutrino temperature, i.e, T^, and 0e are given by 



^ \9tA 9^ 



9'^Ta 



1 + cosh (6* - 
sinh(6' — 0e) 



1 + cosh(6' 
sinh(^ + 



dT, 9T, 



9 fniv) ^ 2^41 
06. ^ 9 



(1 + COsh(^ - (f)e)f (1 + COsh(e + 0e))^ 



1 



1 + cosh(6' 



1 + cosh(6' 



sinh(6' — 0e) _|_ sinh(6' + 0e) 



(A14) 



(A15) 



(A16) 



[1 + cosh{9 - (pe)f (1 + cosh{9 + (pe))\ 
where we used d^y/dT^ = 2y/Tj,. 

The second variable is a perturbation in the total energy density of electron and positron 
induced by i? 7^ 0, i.e.. 



A(pe- + Pe+] 



127r2 



dvfpiv), 



(A17) 
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where 

fpiv) = 7' 



+ 



1 



2(1 + cosh(^ + (f)e)) 2(1 + cosh(^ - 0e)) ^ V 1 + e^+'f'- 1 + e' 



1 1 
+ 



(A18) 



was defined. Partial derivatives of tliis function with respective to Tg, T^, and 0e are 
given by 



dTo en 



sinh(^ + 



sinh(0 — 



1 



2(1 + cosh(^ + 2(1 + cosh(e - 0,))^ 0^ \1 + e''+'^- 1 + e' 



1 



:i + e^+<^'= 



:i+e« 



1 



+ 



1 



2(1 + cosh(^ + 2(1 + cosh(^ - 0^)) ^ V 1 + e^+'^^ 1 + e^ 



+ 



1 



(A20) 



sinh(6' + 



sinh(6' 



2(1+ COsh(^ + (f)e)f 2 (1 + COsh(^ - (f)e)f 



+ 



:i+e^+^^: 



:i+e« 



(A21) 



The third variable is a perturbation in the total pressure of electron and positron, i.e., 

4 



A(Pe-+Pe 



dvfpiv)^ 



(A22) 



where 



/P(^) - 7M ^ 



1 



+ 



1 



2(1 + cosh(e + 0e)) 2 (1 + cosh(^ - 0e)) 



_2--rf/9^ 

was defined. 
3. density-temperature relation 



1 1 

+ 



1 + e^+<^- 1 + 



(A23) 



In the BBN code 



56|, derivatives of 0e with respect to Tg, r = r(T,^) = log(a^) with 



a the scale factor of the universe, and S = J2i ZiYi with Zi the charge and Yi the 
number ratio of nuclide i to total baryon, respectively, are calculated and used. In the 
calculation, we use the following equation for charge conservation in the universe: 



Y 



he 



^=^K-(5)-ne+(5)] 



TT 



(A24) 
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where ng^{B) = neT(O) + AneT{B) is the number density of in an environment of 
magnetic field B, Na = 6.02 x 10^'^ is the Avogadro's number, and k is the Boltzmann's 
constant. 

The left and right-hand sides are denoted as = N{Tg, T^^ipe) and M = M{Tg, r, S), 
respectively. Taking derivatives of both sides with respect t o Tg , r and S, three partial 



derivatives are obtained as in the case of no magnetic field 



dM 



dM 



dr 



Tg,S 



ON ON. 
+ 



r,S 



dN dT^ ^ dN 



dTj, dr d(j)f, dr 



dM 



dS 



dN . 



dr. 



dr 



dS 




dN 



(A25) 



+ 



Tg,S 



dTg 

dN 

(A26) 
(A27) 



Tg,r 



Tg,r dS 

where we used dT^/dr = —T,y/3. The above derivatives are estimated utilizing Eqs. 
(IAI2MAT6I) . and used in estimation of time evolution of the chemical potential param- 
eter (hp. 



dpe/dTg is given by 



dpe dpe , dpe 



dTg dTg 

where dTy/dTg can be described as 

dT, 



dpe dTy 

d(j)e dTg dTy dTg ' 



+ 



dTg 



Ty dr 
3 dTg 



(A28) 



(A29) 



The conservation of energy for mixed matter of 7, i^'s, and baryons leads 56|| to 

(A30) 



dr 


d-P-y _|_ 

dTg 


f dpe _|_ dpe dtf>e ^ 

\dTg ^ d(j)e dTg) 


+ 


dph 
dTg 






dTg 


Pj + ^+pe + ^ + 


f b 1 1 ( dp^ 
IF dr/dt \ dt 


Tg 


_|_ dpe 
' dt 


.) 


T„ dpe 

3 



where p^ and pb are energy densities of photon and baryons, respectively, and P-y and 

Pb are pressures of photon and baryons, respectively. 

Time derivative of the electron and positron energy density is given by 



dpe 



dpe 



Tg 



dr I dt dt 

Using Eqs. (lASOp and f lA31D . we obtain 
dr 



+ 



dS 1 



dpe 

dr ' dS dt dr/dt I ' dT, 



+ 



T 

T 



(A31) 



dTg 



dpi 1 
dTg ^ 


' dpe _|_ 8pe d<j>e ^ 
[dTg d4,e dTg) 


_|_ c'Pb 
dTg 


n ^ P-l ^ n ^ Pe ^ Ph ^ 1 ^Pb 
P7 ^ c2 ^ Pe c2 ^ c2 ^ dr/dt dt 


1 dpe_ 
Tg 9<j>e 


(d4>e 1 dtjledS 1 ^ 

\ dr ' dS dt dr/dt ) 


2rTi dpe 

Z-^^dT^ 
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(A32) 



The derivatives of pe [cf. Eqs. f lA17tiA2ip ] and pressure value Pe [Eq. f lA22tiA23p ] are 
input in this equation. 

It is clear that a PMF enhances energy density of e^+e~ [Eqs. flA17|) and flAlSp ]. The 
PMF, however, does not work on e"*^, so that the enhanced energy of realizes by 
an energy transfer through interaction with thermal bath. When the temperature de- 
creases, the energy gain of decreases because of the weakening of PMF. Resultingly, 
this loss of energy gain heats the thermal bath. This effect is taken into account in Eq. 
f R32|l |86|. 

4. initial value of electron chemical potential 

Using Eqs. f lA12p . flA24p and the fact oi S = in hot environments before the 
nucleosynthesis, we obtain an equation for initial value of electron chemical potential, 
i.e., 



sinh 9 



■^ -1 



cosh I 



(A33) 



where L{z) = K2{z)/z 56|| was defined with K2{z) the modified Bessel functions. 
5. note 

Some transformations in equations are used in order to avoid appearances of divergences 
in a numerical calculation. They include 1/[1 + cosh(^ — (pe)] — + cosh(^ + 0e)] = 



[2sinh0e/(l + cosh(e - (p^Ml - e-2^)/(e'^'= + 2e 



2(1 



)/[e^-(l + e"^') ] and e7(l + e^^ = l/[e7l + e-^)^] 



e sinh x/(l + cosh x)^ 

-x\2-\ 



Appendix B: Effect on nuclear reaction rates 



It was claimed that relatively weak magnetic fields reduce collision rates of nuclear re- 
actions by a factor of two through alignments of directions in which charged nuclei can 



move 



87 



The claim is, however, wrong [89|, |90[ since it was based on a completely- 
mistaken treatment of 2D space perpendicular to the field direction, and cases of i? = and 
B are not treated consistently. 
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Precise reaction rates are derived below. The point is that a target is hit by projectiles 
coming from all directions although the distribution function of projectile charged particles 
are quantized. 



1. Nuclear distribution function 



Magnetic fields would affect nuclear reaction rates through a discretization of momentum 
on the plane perpendicular to the field direction (z-axis). The Zeeman splitting also realizes 
in magnetic fields. It is, however, neglected in this study since situations of large fields are 
eventually excluded from light element abundances deviated through their effects on the 
cosmic expansion rate. Energy levels of charged nuclides are then given by 



Pz 



ZeB(2n + l) 



1/2 



m H 

2m 



pI + ZeB{2n + 1^ 



(Bl) 



where m and Z are the mass and the charge number of a nuclide, respectively. 

The nuclear distribution function is discretized similarly to the case of electron as 



9 



-d?p e" 



oo „ 

ZeBY. J 



9 



-dpz e 



-{E„-fj.c)/T 



where g is the statistical weight, and fiQ is the chemical potential of the nuclide. The number 
density of nuclide i in a Landau level nu and in a momentum range between p^i and Pzi + dpzi 
is given by 



ni{nu,Pzi)dp z 



(2vr) 



'—ZicBexp 



fJ'Ci -m- ZieB{2nu + 1)/ (2mi 



T 



exp 



Pz 



2m,T , 



dpzi- 
(B3) 



The total number density is given by 



ZicB exp 



(2vr) 

l/2Tnl/2 



fJ'Ci -m- ZicB/ {2mi 



T 



exp 



Pz 



2mjT 



dpz 



23/2^^3/2 



Z-fiB exp 



T 



where 



Ui 



rriiT 



1 — exp 



ZicB 

TTLiT 



(B4) 
(B5) 



was defined. 

The fraction in number density is derived from Eqs. flB4p and fIBSP as 

ni{nu,Pzi)dpzi 1 _i / Pli nuZicB' 



21/2^1/2 1/22-1/2 



exp 



27X1 jT 



rriiT 



dpz 



(B6) 
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2. Rates of reactions between charged particles 



The thermal average of reaction rate is described as 



{av) 



oo oo „ 

— XI H / ni{nLi,Pzi)dpzi n2{ni^2,Pz2)dpz2 (y{E)v 



dcos6 



nin2 



(B7) 



"L1=0 riL2=0 



where v is the relative velocity of nuclides 1 and 2, E is the kinetic energy in the center of 
mass (CM) system, and 6 is the angle between momentum vectors of nuclides 1 and 2 on 
the plane perpendicular to the magnetic field. 

The velocity vectors of nuclides i is described as Vi = {v^i, v^i). The angle is then given 
by cos 6* = v^i ■ V ^2 / {t^ ±iv ±2) with amplitudes of the vectors v^i = \v^i\. The relative 
velocity is given by v = [v'j_i — 2v±iv±2fx + v'j_2 + {vzi — ^^2)^]^^^ with /i = cos6'. The CM 
kinetic energy is -E = with the reduced mass /Xrcd = mi'm2/ {rrii + 7712). 

Substituting Eq. (IB6P in Eq. (IB7p . we obtain an equation, 



00 00 



{av) 



20 1/2 l/2r7-i 

"^vrm^ m2 I 



-Ui ^U2 ^ 



J2 exp 

nLi=0 nL2=0 



T 



nil 



m2 



X 



/oo pco rl i 

dpzi / dpz2 / c^/i exp - 



2T 



a{E)v. 



(B8) 



Momentum variables, i.e., Pzi and p^2; are transformed to the CM momentum and the 
relative momentum. Integration over the CM momentum is performed in the equation, and 
we obtain. 



{av) 



1/2 



21/27^1/2 7^1/2 

00 00 
nLi=0nL2=0 



1 — exp ( I 



1 — exp 



niiT 
mi m2 J T 



ZoeB 



m2T 

J dvzr J dfi exp 



2T 



a{E)v, 
(B9) 



where v^r = Vzi — Vz2 is the relative velocity in the direction of the field. 

Velocities on the plane perpendicular to the field are discretized as = {2nu + 
l)Z.ieB /mf. The relative velocity can then be described by 



eB 



(2nLi + l)Zi , (2nL2 + 1)^2 _ ^^ V(2^li + l)^i (2nL2 + 1)^2 



ml 



mim2 



1/2 



(BIO) 
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3. Rates of reactions between a neutron and charged particles 



In reactions of neutron, discrete momenta of only charged particles are taken into account. 
The rate is described as 

(-) = t^^ /n.(nL.P.)*. / {j^f i'p exp (-^) .(B)., (BU) 

Variable transformations from Pzi and Pzn to PzG (CM momentum) and Pzv (relative 
momentum) are performed, and an integration over p^G is computed. The reaction rate is 
then rewritten to be 



{av) 



1 



-1 ^ exp 



nLi=0 



ni^iZicB 



/oo r\ rco I 

dpzv dfi dp^n dpyn exp - 
-oo J —1 J—oo \ 



rriiT 



Pzr Pxn ~^ Pyn 



2/ii 2m„i J 

We perform an integration over azimuth angle on the {pxn,Pyn) plane and trivial transfor- 
mations from momenta to velocities, and obtain an expression, i.e.. 



(av) 



1/2 
A'rcd 



21/27^1/2 2^3/2 



1 — exp 



ZieB 
rriiT 



nLi=0 



roc roc rl I 

Jo Jo J 1 ^'^ i ' 



niiT 



2T 



where the relative velocity is given by 



a{E)v, (B13) 



(B14) 



The kinetic energy in the CM system is given by = 

The sum in the reaction rate is transformed to an integration using the Euler-McLaurin 
formula. The integration form of reaction rate is given by 

1 



{av) 



1/2 



21/2^^1/2 7^3/2 



1 ( ^1^^^ 



x< / dni^i exp (- ^^^^]^^ ) [ dv^rexp 
10 \ mil J Jo 



2T 



^ dv^n V±n exp ^- 



m„.v 



2T 



d^ a{E)v 



dv±n v±n exp ^- 



mr,.v 



2T 



dfi [a{E)v 



nLi=0 
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1 ZieB 



dvzv exp 



12 m\ Jo 

roc 

X / dv±n V±n exp 

Jo 



/w4 \ 

2T ; 

2T 



/ dfi fi [1 —fi 

J~i \ v^i 



d[a{E)v] 



(B15) 



riLl=0, 



Two terms scaling as eB are induced by magnetic field 5, and they disappear in the limit 
of no field, i.e., 5 = 0. 



a. Be(n,p)'^ Li 

We check the reaction rate of ^Be(n,p)^Li for example. Rates are calculated with Eq. 
(lB13p . The masses of neutron and ''Be nucleus is m„ = 0.939565 GeV and mrBc = 6.534184 



GeV 



9l|. Although the cross section would be changed in magnetic fields through the 



momentum quantization, that effect is neglected and cross section values in no fields are 



taken from Ref. [6l[ approximately. 

Figure m shows calculated rates of ^Be(n,p)^Li as a function of Tg. Solid dark lines corre- 
spond to cases in which magnetic fields decrease with time because of the cosmic expansion 
as -B oc T^. Field amplitudes are B/{BqT^) = 10^, 10"^ and 10^, in order of increasing line 
width (from the top to the bottom), respectively. Dashed lines, on the other hand, corre- 
spond to cases of constant magnetic field of B/Bc = 10^, 10'^ and 10^, in order of increasing 
line width (from the top to the bottom), respectively. Solid gray line corresponds to the 
rate in no magnetic field. The uppermost solid and dashed lines are hardly distinguishable 
from the solid gray line. 

In the attenuating magnetic field case (solid dark lines), the field effect is larger in higher 
temperatures. The discretization effect is roughly determined from the index factor, i.e., 
ZieB/{miT), in the exponential in Eq. flB13|) . The field amplitude (oc T^) decreases faster 
than the temperature in the universe (oc T^) does. The index is, therefore, larger in higher 
temperature. In the constant magnetic field case (dashed lines), the field effect is larger in 
lower temperatures conversely. 

In somewhat high magnetic field, the minimum energy of nuclear Landau level is higher 
than the thermal energy of the universe, ~ T. The average of rate then receives a contribu- 
tion from large CM energies which originate from large relative velocities [Eq. (]B14p ]. Since 
the reaction rate at higher energies is roughly smaller as for the reaction ''Be(n,p)^Li, the 

21 



0.1 



10 



FIG. 4. Reaction rates of ^Be(n,p)^Li as a function of the temperature Tg = T/(10^ K). Solid 
dark lines correspond to cases of attenuating magnetic field for B / {BqT^) = 10^, 10"^ and 10^, in 
order of increasing line width, respectively. Dashed lines correspond to cases of constant magnetic 
field for B/Bq = 10^, 10'^ and 10^, in order of increasing line width, respectively. Solid gray line 
corresponds to the rate in no magnetic field. 

existence of field decreases the reaction rate. The lowest solid dark line in Fig. H] has two 
bumps at 1 ^ Tg ^ 2 and 4 ^ Tg ^ 7. These bumps correspond to peaks in reaction rates 
produced by the S"*" resonant states of ^Be* at resonance energies Ej- = 0.33 MeV and 2.66 
MeV {92!. At Tg = 1.75 and Tg = 5.2, the minimum CM energies of the ground Landau 




level, i.e., E^^^ = jj^^^ZieB /{2ml) [cf- Eq. flBTi]) ]. are 0.3 MeV and 2.7 MeV, respectively. 

As observed above, magnetic fields can affect thermonuclear reaction rates. Large am- 
plitudes of magnetic fields as assumed in Fig. H] are, however, excluded from incredibly fast 
expansion of universe (see Fig. [2]). The effect of magnetic field on nuclear reaction rates can 
thus be neglected. 
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